Determining gap nodal structures in Fe-based superconductors: angle-dependence of the low 
temperature specific heat in an applied magnetic field 
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Since the discovery of high-Tc LaOi-a;Fa;FeAs and other such systems based on FeAs layers, several propos- 
als have been made for the superconducting order parameter Ak, on both phenomenological and microscopic 
grounds. Here we discuss how the symmetry of Ak in the bulk can be determined, assuming that single crystals 
will soon be available. We suggest that a measurement of the dependence of the low temperature specific heat on 
the angle of a magnetic field in the FeAs plane is the simplest such method, and calculate representative specific 
heat vs. field angle oscillations for the various candidate states, using a phenomenological band structure fitted 
to the DFT Fermi surface. 



PACS numbers: 74.70.-b,74.25.Ha,74.25.Jb,74.25.Kc 

The recent discovery of superconductivity with onset tem- 
perature of 26K in LaOi_2;Fa;FeAs [1] was followed rapidly 
by the development of materials with Tc up to ~50K[2-10], 
which possess a similar structure but where La has been re- 
placed by Sm, Pr, Nd or Ce. Common to all such materi- 
als is an electronically layered structure, where according to 
electronic structure theories a rare-earth oxide layer where F 
substitutions for O dope an FeAs layer. The iron atoms are ar- 
ranged in a simple square lattice, separated by arsenic atoms 
above and below this plane. The FeAs complex provides the 
bands at the Fermi level, and the Fermi surface consists of 
sheets around the F point and the AI point of the Brillouin 
zone[l 1-16]. 

Beyond this general initial consensus on the commonalities 
of the different materials, electronic structure calculations dif- 
fer on the details of the ground state and the band structure 
near the Fermi level. Both paramagnetic and antiferromag- 
netic (sublattice and linear SDW) ground states have been re- 
ported, with some authors claiming that the system is close to 
a Mott transition and also possibly to a ferromagnetic state. 
Crude support for the proximity of competing magnetic states 
is provided by the known helimagnetism in the layered iron 
monoarsenide system. 

Within weeks of the discovery of the LaOi_2;F3;FeAs sys- 
tems, theoretical analyses of various possibilities for the 
mechanism of superconductivity and the symmetry of the 
superconducting order parameter have appeared[12, 13, 17- 
19]. Eliashberg style calculations based on density func- 
tional theory (DFT) determination of electron-phonon cou- 
pling constants [20] suggest that conventional electron-phonon 
interactions are not sufficient to generate the observed transi- 
tion temperatures. Thus several authors have discussed elec- 
tronic pairing mechanisms of the spin fluctuation type[12, 13, 
17-19], but disagree about the symmetry of the ground state, 
apparently because of the details of the electronic structure 
used as an input to the calculation. Given the past history of 
theoretical approaches to unconventional superconductors, it 
may be some time before a consensus on the coiTect micro- 



scopic approach is forged. 

In the intervening period, it would clearly be useful to have 
some information on the symmetry of the order parameter 
to guide such theoretical discussions. Evidence for nodes in 
the order parameter has already been provided by point con- 
tact tunnelling[21], which has reported a zero bias state in a 
series of relatively high-transparency junctions, and specific 
heat measurements in a magnetic field H[22], which indicate 
a Cv/T ^ ^/H term similar to that predicted by Volovik 
for a d-wave (or, more generally, nodal) superconductor. Be- 
cause the current experiments have been performed on pow- 
dered samples, however, the distribution on the Fermi surface 
of order parameter nodes, which could provide some infor- 
mation on the symmetry of the pair state is not yet deter- 
mined. In addition, the point contact measurements probe 
only the superconducting state at the surface, whereas ide- 
ally one would prefer to extract information on the bulk su- 
perconducting state. When single crystals are produced, it 
will be possible to perform what is possibly the simplest bulk 
probe of the distribution of gap nodes, a measurement of the 
specific heat of a sample in the presence of a field in the 
FeAs plane as a function of its angle relative to the crystal 
axes. The superflow field in the vortex state of the type-II 
superconductor is known to "Doppler shift" the energies of 
quasiparticles, changing their local occupation and giving rise 
to a residual density of states[23] which depends on the an- 
gle the field makes with the nodes[24, 25]. This method of 
"nodal mapping" was proposed in the context of the cuprates, 
where the experiment is difficult due to the large phonon back- 
ground, but has found more fruitful application in lower- Tc 
materials [26]. 

In this paper we calculate the specific heat oscillations with 
magnetic field angle to be expected in the presence of a va- 
riety of candidate superconducting pair states. Rather than 
tie ourselves to any particular microscopic electronic structure 
calculation, we use a phenomenological two-band model [27] 
which captures the essential qualitative features of the bands 
near the Fermi surface. We find that various extended- s like 
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FIG. 1: (Color online) The two hybridized bands of our model. The 
white region depicts the narrow energy window around the Fermi en- 
ergy (dashed line) where the model reproduces the semiquantitative 
aspects of the LDA band structure. The momenta refer to points in 
the 2D effective large Brillouin zone (see Fig. 2). 




states can be distinguished from, e.g. rf-wave or p-wave like 
states by the positions of their nodes. There are also cases, 
however, where nodes lie in positions on the Fermi surface 
where k„ and the Fermi velocity vp are not parallel. In this 
case the minimum of the specific heat does not correspond 
precisely to the nodal position and the structure of the set of 
minima must be examined in detail. 

Effective band structure. The crystal structure of LaOFeAs 
consists of alternating layers of FeAs and LaO. Density func- 
tional theory (DFT) calculations show that the energy bands 
crossing the Fermi level can be assigned to the Fe 3d and the 
As 4p-orbitals [11-17]. Thus, to describe superconducting 
properties we can consider the LaO layers mainly as spacing 
layers and possible charge reservoirs. The FeAs layers can be 
further subdivided into a square Fe lattice with an Fe-Fe spac- 
ing dpe-Fe = 2.82 A and an As square lattice displaced by a 
vector (1/2, 1/2) in the x-y-plane to a position in the center 
of the Fe squares. Additionally, the As atoms are displaced 
alternately above or below the Fe plane leading to a pyrami- 
dal Fe-As configuration. Due to the alternating sign of the As 
^-displacements the primitive unit cell of LaOFeAs contains 
two Fe and two As atoms. The axes of the corresponding Bril- 
louin zone (BZ) are aligned in the next nearest neighbor Fe-Fe 
direction and the BZ has a size of 27r/a x 27r/a. However, due 
to the high degeneracy of the two As positions we will treat an 
effective model consisting of a smaller unit cell having only 
one Fe and one As atom. This leads to a larger effective BZ 
that has axes that are aligned to the nearest neighbor Fe-Fe 
direction. In this case the real BZ occupies a diamond shaped 
region within the effective BZ. 

Since the thermodynamic properties of the superconducting 
state are governed by low lying quasiparticle excitations, we 
restrict our considerations in the following to the region where 
the energy bands cross the Fermi level, forming the different 
Fermi surface sheets in the BZ. To simplify the rather compli- 
cated band structure, we use a two band model that takes only 
the iron d^z and dyz orbitals into account [27]. Here the basic 
symmetry of the hopping parameters is determined from the 



FIG. 2: (Color online) The different FS sheets in the large effective 
BZ calculated within our two-band model. The two hybridized bands 
result in two sets of Fermi surface sheets, centered around the F point 
(blue) and the AI point (red) of the real BZ (dashed black line). Back- 
folding of the large into the small BZ produces the dashed sheets in 
the small zone. The black crosses show the FS of the paramagnetic 
ground state determined by DFT [14]. 



direct overlap of the Fe d orbitals as well as from the hopping 
mediated by the As p orbitals. The model neglects contribu- 
tions from other orbitals, e.g. hybridization due to the other 
Fe d orbitals, and the hopping parameters are adjusted to give 
the generic form of the Fermi surface sheets determined by 
bandstructure calculations. 

Calculating the hopping from the direct overlap of the Fe 
dxz and dyz orbitals as well as the effective hopping in second 
order perturbation theory on the path Fe-As-Fe, taking the As 
Px, Py and pz into account, leads to a tight binding Hamilto- 
nian with nearest and next-nearest neighbor hopping between 
the same orbitals and a next-nearest neighbor exchange hop- 
ping between the two bands. Due to the choice of the orbitals 
there are different nearest neighbor hopping values <i and t2 
for hopping in the x and y directions in one band which are 
interchanged in the other band. The intraband next nearest 
neighbor hopping is the same for both bands and both direc- 
tions, while the interband hopping t4 has a different sign for 
the (1,1) compared to the (1,-1) direction. After the usual 
Fourier transformation we can write the intraband energies in 
momentum space as 

£11 — — 2t I COS kx — 2t2 COS ky — At^ COS kx COS ky (1) 
£22 — ~2t2 COS kx — 2ti COS ky — Ats COS kx COS ky (2) 

and the interband exchange energy is 

£12 = £21 = -4*4 sin kx sin ky (3) 

Taking the hybridization of the two bands into account one 
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FIG. 3: (Color online) Candidate order parameter states considered 
in this work, a) d^2_y2 state; b) dxy state; c) extended-s state from 
Ref. [13] ; d) generalized extended-s; e) state (Ak oc sin k^)', f) 
Px + ipy state (Ak = sin k^+i sin ky) [12]. The dashed orange line 
denotes the small Brillouin zone and the green line denotes locus of 
gap nodes. 



finds for the two bands 
1 



£11 + £22 - V (£11 - £22)^ + 4e^; 
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£11 + £22 + V (£11 - £22)^ + 4ef; 



(4) 
(5) 



The hopping parameters ti and the chemical potential Ep can 
be used to fit the FS of the paramagnetic ground state found 
within the DFT calculations [14]. We find a reasonable agree- 
ment of the FS for the following values ti ~ —1.2, ^2 ~ 1.35, 

= —0.8, ^4 — —0.8 and Ep — 2. These values lead to the 
bandstructure shown in Fig. 1. 

Order parameters. Because of the small coherence length 
of 0(30 A), these systems are strongly type-II and it may 
be appropriate to think of the range of the pairing interac- 
tion as being very short, of order the lattice spacing. Gen- 
erally speaking, order parameters involving pairing on nearest 
neighbor sites are also those which have been proposed for 
these systems. We therefore consider as representative can- 
didates the states listed in Fig. 3, proposed by various au- 
thors, beginning with nearest-neighbor a) dj.2_y2 state (Ak oc 



cos kx — cos ky)[19] ; b) d^y state (Ak oc sin sin ky); and 
c) extended s-wave state (Ak oc cosfcx + cosfcj^) [13], and 
e) a Px-wave state. The extended s-wave state shown in c) 
changes sign on the Fermi surface of the model system, as 
seen in the Figure. On the other hand, its nodes are located at 
the 45° directions relative to the sheet center at the IVI point, 
which is not generic for a state with s (Aig) symmetry. We 
therefore show in Fig. 3d) the result of adding to this state a 
higher order s-harmonic Ak oc (1 — 'y){coskx + cos ky) + 
7(cos2fca: + cos2A;j,) with 7 = 0.05. For example, the RPA 
spin fluctuation calculations of Kuroki et al. appear to lead 
to a more general extended s-wave state. Similarly, it can be 
seen in Fig. 3d) that the points where the nodes cross the 
Fermi surface sheets are away from the 45° directions (rela- 
tive to the center of the sheet on the zone face). Finally, we 
show in Fig. 3f) the nodeless px + ipy state proposed by Xu 
etal. [12]. 

Specific heat. To get a qualitative understanding of the spe- 
cific heat oscillations as a function of the rotation angle of an 
in-plane magnetic field at low temperature it is sufficient to 
study the spectrum of low energy excitations. To calculate the 
spectrum in the vortex state we want to follow a semiclassical 
approach, that neglects the core states and considers only the 
shift of the quasiparticle energies of the extended nodal states 
in the presence of a magnetic field. Following [25] we approx- 
imate the vortex lattice using a circular unit cell with radius R 
and winding angle (3. Then the Doppler shifted quasiparticle 
energy is 



oE^ ' = m-vy-Vs — 



P y COS a — Vp^ sm a 1 sm p 



(6) 

Here v^^ denotes the Fermi velocity on band i, is the gauge 
invariant expression of the quasiparticle flow field around the 
vortex core and a is the angle between the magnetic field and 
the a;-axis of our coordinate system. The dimensionless radial 

variable p ~ r/ R and v^Pxjy components of the Fermi 

u\ 

velocity calculated from Ve). , normalized by a Fermi surface 
averaged value of Vp . The energy scale En associated with 
the Doppler shift is 



E 



(7) 



where a is a geometric constant characteristic of the vor- 
tex lattice, $0 is the flux quantum, and Vp is an averaged 
Fermi velocity on band i determined from DFT calculations. 
This procedure of normalizing the Fermi velocities prevents 
us from overestimating the differences in the energy gradients 
at the Fermi level of the simplified two-band model. Using 
the Doppler shifted energy in a BCS-like density of states we 
can calculate the low energy spectrum as 
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FIG. 4: (Color online) Residual density of states N{uj — 0, H) vs. 
a, the angle H makes with the a;-axis for the different superconduct- 
ing states shown in Fig. 3. N{lo — 0, H) is proportional to the 
linear specific heat coefficient at low temperature. Note all curves 
have been offset by a constant amount for clarity. 



where the angular brackets denote an average over the vortex 
cell (H) and over the Fermi surface (FS), respectively. The 
integral over the vortex cell can be done analytically leading 
to 



i,(i?«/A.)^(4:i 



cos a 



vl,' sm a 



(9) 

The last average is to be performed over the different Fermi 
surface sheets in the unfolded Brillouin zone leading to an 
oscillation of the low energy spectrum as a function of the an- 
gle of the applied magnetic field. These oscillations can be 
directly determined by low temperature thermodynamic mea- 
surements, like specific heat or the thermal conductivity. 

In Fig. 4, we show the residual angle-dependent density of 
states, or linear specific heat coefficient as a function of field 
angle 9. For the most part, one expects fairly straightforward 
generalizations of the results for a circular Fermi surface [25], 
as seen for the and d-symmetry states: minima in the spe- 
cific heat at low temperatures T <^ E\j lie at the expected 
nodal positions. The nodeless + ipy state produces no 
Volovik effect, is therefore not plotted in Fig. 4 and is ap- 
parently not a candidate for the Fe-based materials. In the 
extended-s cases, some interesting points arise. It is seen from 
Fig. 3c) that in the simple extended-s case, the nodes are lo- 



cated along the 45° directions. Nevertheless the minima in 
Fig. 4 are slightly displaced symmetrically with respect to 
these nodes; this is due to the fact that the M sheets are ellipti- 
cal, with the consequence that the Fermi velocities are not par- 
allel to nodal k„ measured from the sheet center When there 
are higher harmonics, such as in the generalized extended s- 
wave case, the nodes themselves actually are displaced from 
the 45° directions. Thus a measurement of this kind can iden- 
tify an extended s state by the displacements of the minima, 
but a direct correspondence with the nodal positions requires 
a precise knowledge of the underlying Fermi surface. 

Conclusions. In this paper, we have proposed that the mea- 
surement of specific heat oscillations as a function of the mag- 
netic field angle in the FeAs plane of the new iron-based su- 
perconductors could be the simplest and most straightforward 
bulk determination of gap symmetry, once single crystals or 
highly oriented powders are available. To simplify the calcu- 
lation, we used an effective two-band model with parameters 
chosen to reproduce the DFT Fermi surface. We then calcu- 
lated the low-temperature linear term in the specific heat to be 
expected as a function of field angle for a variety of candidate 
states. The elliptical Fermi surface pockets near the M points 
introduce some interesting complications in the problem rela- 
tive to the usual picture of C\/(H) oscillations over the field 
angle. 
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